TDR-269fS48i0-22)-2 


W/E  t>4  -  1007 
Do  v  SB  PUB  control  only) 


At)60/JjLb 


NON  EQUI  LI  B  RJ  UM  SODI UM  IONI 7J  FI  ON 
IN  LAMINAR  BOUNDARY  LAYERS 


14  April  1964 


Prepared  by 
John  J  Kane 


AEROSPACE  CORPORA TION 
SAN  BERNARDINO  OPERATION 
San  Bernardino,  California 


Contract  No.  AF  04(695) -269 


Prepared  for 


COMMANDER 

DQU A  R TE RS ,  BA  L  U S  TI C  S YS  TE MS 
AIR  FORCE  SYSTEMS  COMMA N 
UNITED  STATES  AIR  FORCE 


DIVISION 


N crio n  kAir  F o r c e  B a s e 


TDR-269(S4810-22)-2 


(This  abstract  is  unclassified) 
ABSTRACT 


The  injection  of  sodium  as  a  trace  contaminant  into  a  laminar 
boundary  layer  is  considered.  Three  ionizing  reactions  are 
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Comparisons  are  made  with  the  predicted  clean  air  variations 
with  the  same  parameter.  Stream  tube  predictions  (neglecting 
diffusion)  are  also  made  and  compared  with  the  more  exact 
results.  The  relative  ionization  of  clean  air  and  sodium  is 
assessed  as  a  function  of  altitude  and  amount  of  injection 
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I.  INTRODUCTION 

The  effects  of  mass  transfer  in  high-speed  flow  have  been  a  problem  of 
continuing  interest  in  the  design  of  vehicles  subjected  to  intense  aerodynamic 
heating.  To  gain  heat  protection,  gases  or  liquids  can  be  actively  injected 
into  the  boundary  layer  at  controlled  rates,  a  method  known  as  transpiration 
cooling  or  mass  transfer  cooling.  Foreign  species  can  also  be  injected  into 
the  boundary  layer  at  a  rate  not  as  positively  controlled,  if  the  vehicle  is 
provided  with  an  ablating  surface. 

The  injection  of  species  that  do  not  chemically  react  with  the  air  in  the 

1  2  ctnd>  3 

boundary  layer  has  been  studied  in  some  detail  ’  ’  .  However, 

because  the  boundary  layer  equations  become  much  more  complicated  and 
because  of  a  lack  of  sufficient  chemical  kinetic  data,  little  work  has  been 
done  on  mass  transfer  into  a  chemically  reacting  non -equilibrium  boundary 
layer.  One  problem  of  this  type  which  is  quite  significant  a"  \  which  lends 
itself  to  a  ready  solution  is  the  problem  of  trace  contaminants. 

An  ablation  heat  shield  can  in  manufacturing  or  handling  be  contaminated 
with  very  small  amounts  of  certain  compounds.  While  insigr ' 'icant  with 
respect  to  their  total  weight,  they  can  be  significant,  if  they  are  highly 
ionizable,  and  when  injected  in.o  a  boundary  layer  become  dominant  in  the 
production  of  electrons.  The  observable  characteristics  of  a  vehicle  can  be 
determined  by  these  trace  contaminants. 

This  report  evaluates  the  effect  of  trace  contaminants  on  the  boundary  layer 
with  mass  addition.  It  assumes  that  the  trace  contaminants  do  not  affect  the 
velocity  and  enthalpy  profiles  and  these  profiles  are  determined  by  the  gross 
amount  of  blowing  into  the  boundary  layer.  Similar  solutions  are  then  derived 
for  the  distribution  of  ionized  trace  contaminants  through  the  boundary  layer. 
Several  types  of  ionizing  reactions  are  considered.  It  is  assumed  that  the 
degree  of  ionization  is  much  Jess  than  equilibrium  and  is  determined  by  the 
reactions  and  rates  of  Bortner. 
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BASIC  EQUATIONS  AND  SOLUTION 


Consider  a  contaminant  being  injected  into  the  boundary  layer  along  a  cone 
in  small  enough  amounts  so  that  it  does  not  significantly  alter  the  velocity 
and  enthalpy  profiles.  Then  one  need  only  consider  the  species  equation. 

Consider  the  ions,  atoms,  and  electrons  associated  with  the  trace  contaminants 
as  one  species.  The  governing  equation  for  the  distribution  of  this  species 
through  the  boundary  layer  is  then: 

d  f  1  dC  .  .  dC  _  n  ,,, 

dtj  ^  S?  drf  j'  +  £  dT  '  °  (1) 

where  f(’l)  represents  a  known  velocity  function  determined  from  the  gross 
blowing  into  the  boundary  layer^. 

Boundary  conditions: 

At  r\  =  0,  C  =  C 

w 

At  rj— ►  oo ,  C  “  0 

i  cw  _  c 

Let  be  constant  and  S  =  — ^ -  (2) 


Then  Equation  (1)  becomes 

l  d2S  .  ,  dS 

Sc  ,2  *  dq  " 

d  q 

Boundary  conditions: 

At  t)  s  0,  S  =  0 

At  t]— oo,  S  =  1 
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From  the  momentum  equation  for  constant  f  it  is  readily  seen  that: 

Sc 


exp 


L  f 11  f  .  1  (  fn  (0)  I 

i  Sc  ( l  T  \ 


(5’ 


Substituting  into  Equation  (4),  we  get  a  more  convenient  form  of  solution 

S  =  -—4 - -  (6) 


f  /'!D. 

Sc 
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i 

Sc 

dri 

This  represents  the  distribution  of  total  injected  trace  species  (atoms,  ions, 
and  electrons)  across  the  boundary  layer. 

Consider  now  the  production  of  ions  through  the  boundary  layer  where 
the  source  term  can  be  represented  as  a  function  of  t)  only.  The  governing 
species  equation  is  then 'neglecting non- similar  terms): 


a  d^C  LfdC  ^  2  x  w  n 

s?  7F  - 0 


Boundary  conditions: 

At  T)  S  0,  c  -  C 

At  "n  =  oo,  C  ~  0 

where  w  =  w  (t|) 

The  solution  to  this  equation  is: 
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where 


(9) 


Thus,  if  the  source  term  w  can  be  described  as  a  function  of  q  only,  then 
Equation  (8)  will  determine  the  ion  distribution  through  the  boundary  layer. 

III.  IONIZATION  KINETICS 

Sodium  as  a  trace  contaminant  was  studied.  The  particular  ionization 
reactions  considered  are  listed  below.  Any  other  trace  contaminant  could 
be  studied  with  the  existing  programmed  solution  as  long  as  it  were  ionized 
by  any  or  all  of  these  mechanisms  and  the  corresponding  rate  constants  were 
known. 


The  particular  reactions  included,  which  are  those  considered  most  significant 
for  sodium  ionization^,  are: 


1. 

Na 

+ 

M  -  Na+  + 

e  + 

M 

(10) 

2. 

Na 

+ 

e  -*  Na+  + 

2e~ 

(11) 

3. 

Na 

+ 

0  4  0-  Na+ 

+  e 

+  °2 

(12) 

where  M  refers  to  any  third  body  (except  O,  and  e  when  reactions  2  and  3 
are  included). 

The  forward  direction  is  considered  only  for  each  reaction,  implying  that  the 
reactions  are  not  close  to  equilibrium.  As  a  check  on  the  validity  of  this 
assumption  the  equilibrium  concentration  can  be  computed  for  each  particular 
case  studied  and  compared  with  the  non -equilibrium  prediction  of  this  analysis 
(see  Appendix  A). 
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From  the  described  ionization  reactions  and  the  law  of  mass  action^1,  , ,  an 
expression  for  the  source  term  follows  as: 


w 

P 


where 
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(1  -  C0)  +  k 2  fV  k 
Me„  lc3 


(CQ) 

Mq 


l 


t 

\ 


13, 


p  =  mixture  density 

Mm  =  mean  molecular  weight  (Reaction  1) 

^Na  =  mass  fraction  of  sodium  atoms 
kj  =  rate  constant  associated  with  reaction  j 
Cj  =  mass  fraction  of  species  i 

-  molecular  weight  of  species  i 

To  a  first  approximation,  the  distribution  of  sodium  atoms,  Cj^,  is  given 
by  Equation  (6).  This  will  be  a  very  accurate  representation  unless  the 
degree  of  ionization  is  very  high. 

5 

The  rate  constants  associated  with  each  reaction  are  given  by:  > ,< 


kl 


k2 


3/2  ~ 107»  500  ^  i  i 

=  80  r  '  ^  e  T  fr  (lb-mole)  sec 


11  1/2  ~ 107 , 500  ^ 

1  (10)1A  T1/4  e  T  fr 


=  3. 1  (10) 


(lb-mole)  ^  sec  * 


=  10.3  (10)^  ft''  (lb-mole)  ^  sec  * 


(14) 

(15) 

(16) 


The  mass  fraction  of  electrons  acting  to  ionize  sodium  by  reaction  2  are 
considered  to  be  those  electrons  produced  by  reactions  1  and  3  as  well  as 
the  electrons  associated  with  the  ionized  NO  in  the  boundary  layer.  The 
clean  air  distributions  of  atomic  oxygen  and  ionized  NO  can  be  determined 
by  the  procedure  described  in  Reference  7.  Figure  1  describes  the 
distributions  used  in  this  report. 
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The  source  term  is  thus  completely  described  as  a  function  of  q  except  for  the 
temperature  dependency  of  the  rate  constants. 


For  a  Lewis  number  of  unity  and  a  constant  Prandtl  number,  the  energy 
equation  is  (assuming  similarity): 


f  g‘  + 


(h~  -h 
v  se  sw' 


h  > 

t  Pr 


!  {■  f"  \  =o 


hs  M  "hsw 
where  g  =  t- - — - 

^*00  iAottr 


hs  M 


=  total  enthalpy 


=  h  (q)  + 


U2 

T~ 


prime  denotes  differentiation  with  respect  to  q. 

Boundary  conditions: 

At  q  is  0:  g  =  0 
At  t|  =oo:  g  =  1 

The  solution  to  Equation  (17)  is: 


g(-n)  =  Ci  f  |  £"  (n) 

J  n 


•  T  L  |f"  M>|  "f0'  H(t|)|£"(>l)  “dndn 


Pr  yr\ 


(18)  - 


where 


Pr/"00  i  _Pr 

C,  =  H  (q)  J  f"  (q) 


dq  t 
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and 


H  (*l) 


i  ue2 

^se  '  ^sw 


f‘  (T|S  £•”  (tl>  +  f»  (11)  f"  (Tl) 


(20) 


Equation  (18)  can  be  used  -.o  describe  the  total  enthalpy  distribution  through 
the  boundary  layer  and  m  conjunction  with  the  known  velocity  profile  gives 
the  static  enthalpy  or  temperature  profile.  Thus,  the  source  term  is 
completely  described  as  a  function  of  q  and  Equation  (8)  can  be  used  to  give 
the  ion  (and  electron)  distribution  through  the  boundary  layer.* 


IV.  RESULTS 


Solutions  to  Equation (8) were  computed  for  each  of  the  three  ionizing  reactions. 
The  Schmidt  number  for  the  diffusion  of  the  sodium  atoms  into  the  boundary 
layer,  as  determined  from  Equation  (6),  was  considered  to  have  the  constant 
value  0.5.  The  Schmidt  number  for  the  diffusion  of  the  ion-electron  pairs 
as  determined  from  Equation{8)  was  taken  as  0.25.  The  results  are  described 
in  Figures  2  and  3  which  represent  the  scaling  of  electron  density  versus  the 
parameter  Px/Ue.  The  electron  density  is  scaled  in  the  form  Ne/Pm  Kwn 
where: 

N  =  electron  density  /electrons /cc) 
e 

P  =  pressure  (atmospheres) 

Kw  =  wall  concentration  of  sodium 
m  =  order  of  initial  reaction 


n  -  1  except  when  the  sodium  is  ionised  by  an  electron  which  itself 
is  a  product  of  the  ionizaLon  of  sodium,  in  which  case  n  is  equal  to  two. 

The  non -equilibrium  clean  air7  and  sodium  equilibrium  ionization  levels 
are  included.  As  seen  in  Appendix  B,  the  equilibrium  ionization  scales  as 

(Ne)eq  Ne 

— Y~Ty‘ -  while  the  clean  air  non-equilibrium  ionization  scales  a  s  . 

p1 1  c  K  t 


*This  solution  and  the  equilibrium  analysis  of  Appendix  A  have  been  programmed 
in  FORTRAN  for  the  IBM  7094. 
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Each  curve  represents  the  peak  electron  concentration  generated  by  the 
particular  reaction  in  the  boundary  layer  with  a  certain  peak  temperature. 

For  different  peak  temperatures,  the  curves  would  be  shifted.  To  indicate 
the  degree  of  temperature  dependency,  two  cases  of  different  peak  temper¬ 
atures  were  studied  as  indicated  in  the  figures.  The  peak  temperature 
actually  refers  to  the  temperature  at  the  location  at  which  the  electron  mass 
fraction  is  a  peak  and  is  not  necessarily  the  peak  temperature  of  +he  boundary 
layer.  Blasius1  type  velocity  profiles  were  used  although  provision  was 
included  in  the  program  to  use  the  velocity  profile  associated  with  any  value 
of  the  blowing  parameter,  f(o).  The  atomic  oxygen  and  distributions 

corresponding  to  the  given  wall  and  free  stream  boundary  conditions  were 
used.  (Figure  1). 

The  most  significant  reaction  in  the  non -equilibrium  ionization  of  sodium  is 

the  ionization  by  collision  with  atomic  oxygen.  For  example,  if  the 

contributions  from  each  curve  are  evaluated  for  a  peak  temperature  of  809  i°R, 

-4 

4  |  )  a  wall  concentration  of  10  ,  a  pressure  of  0.  186  atmospheres,  and  a  value 

of  Px/Ue  of  0.262(10)  the  following  is  obtained  (corresponding  to  Q  =  12°, 
h=  115,000  ft,  Uco  =  22,  000  fps,  T^  -  1000°K,  =  10  ^  and  a  cone  length  of  3  ft). 

Table  1.  Electron  Ion:  zation  Contributions 
Curve  Reaction  Ne  (electrons/cm3) 


1 

Atomic  Oxygen 

2.  2  (10)9 

2 

Electron 

2.7  (10)3 

3 

Clean  Air 

8.  8  (10)8 

4 

Equilibrium 

1.2  (10)12 

5 

Electron 

1  (10)3 

6 

Arbitrary  Third  Body 

7.7  (10)5 

7 

Electron 

0.9 

<  »  ) 
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None  of  the  other  sodium -ionizing  reactions  are  within  even  several  orders 

i* 

of  magnitude  of  that  by  which  the  sodium  is  ionized  by  atomic  oxygen.  This 
result  applied  for  all  the  conditions  investigated. 

Table  1  also  indicates  the  significance  of  non -equilibrium  sodium  ionization 
in  comparison  to  clean  air  ionization  for  the  conditions  considered.  If  the 
wall  concentration  had  been  as  low  as  one  part  per  million,  the  sodium 
ionization  would  be  negligible.  However  a  wall  concentration  as  high  as  10-^ 
can  be  quite  significant. 

If  diffusion  is  neglected,  a  simplified  analysis  can  be  accomplished  to 
determine  non -equilibrium  ionization  levels.  This  has  been  done  in  Appendix 
B  and  the  results  are  shown  in  Figure  4  in  comparison  with  the  results  of 
the  analysis  based  on  the  diffusion  equation.  Included  also  is  a  similar 
comparison  for  the  clean  air  computation.  It  is  seen  that  in  both  cases  the 
streamtube  (no  diffusion)  method  over-predicts  by  roughly  two  orders  of 
magnitude , 


To  determine  specific  regimes  in  which  clean  air  ionization  predominates 
over  sodium  ionization,  the  electron  densities  at  various  altitudes  were 
computed  for  the  conditions  of  the  example  of  Table  1.  The  results  are 
described  in  Figure  5.  It  is  seen  that  clean  air  ionization  predominates 
at  higher  altitudes  and  the  sodium  ionization  at  low  altitudes.  The  cross-over 
point,  if  any,  depends  on  the  wall  mass  fraction  of  sodium.  This  change  in 
predominant  ionization  is  due  to  the  fact  that  the  clean  air  ionization  is  a 
first  order  reaction  while  the  principal  sodium  ionization  reaction  is  second 


order. 
5 quar 6 


Thus  as  altitude  decreases,  the  sodium  ionization  increases  with  the 
of  the  pressure  as  well  as  with  (Px/Ue)“  while  the  clean  air  ionization 


is  proportional  to  the  pressure  itself  as  well  as  to  (Px/Ue)^ 


The  above  example,  and  almost  all  the  results  of  this  analysis  are  highly 
dependent  on  the  wall  mass  fraction  of  the  contaminant.  If  this  has  a  value 
of  10  ^  or  less,  the  sodium  ionization  will  be  negligible  to  altitudes  below 
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90.000  feet.  If  it  has  a  value  of  10  or  greater,  the  sodium  ionization  will 
predominate  at  altitudes  below  approximately  200,000  feet.  At  intermediate 
values  the  particular  conditions  will  lave  to  be  investigated. 

The  wall  mass  fraction  of  sodium  c.’n  be  related  to  the  blowing  parameter 
f(o)  and  the  level  of  heat  shield  contamination  by  (see  Appendix  C' 


_  Q  f(o) 

w  -  f(ay-,'£n<oT 


where 


(21) 


Q  =  fractional  sodium  contamination  of  heat  shield  material  (parts 
per  million/ 10^). 


For  a  cone. 


f(o) 


(P  v>  W  V  X 
"\  / Pe  He 


(22) 


If  a  similar  solution  is  to  apply,  then  (pv)  w  /x  r  constant  (23) 

The  results  presented  are  based  on  a  Blasius  profile]"  f(o)  =  oj  ,  implying 
that  very  little  or  negligible  blowing  is  taking  place.  The  Blasius  profile  is 
an  accurate  representation  to  within  approximately  10  percent  of  the  profiles 
for  blowing  up  to  f(o)  =  -0.  10.  Thus  ’he  res  JUs  obtained  using  the  Blasius 
profile  should  be  valid  for  values  of  K  .  Q  up  to  0.2,  corresponding  to  a 
twenty  percent  mass  fraction  of  total  ablauon  products  at  the  wall. 


For  a  case  without  negligible  blowing,  f(o)  would  be  known  as  a  function  of 
altitude.  Equation  (21)  then  would  describe  Kw  as  a  function  of  altitude  and 

o  r*Ar»-A<*nii1iVv7*inm  analtrctc  r  nnl  rl  ho  nprf  nrm  P  H  11  ciyio  thp  DTHTlPr  K  JlTlH 

w-w  - - - - - o  —  ~  r - r —  yfif - 

velocity  profile  for  each  altitude.  ^ 


»  \  1 
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V.  CONCLUSIONS 

The  following  conclusions  apply  for  slender  cones  using  Bortner's  rates 
for  the  ionization  of  sodium. 

1.  Sodium  non- equilibrium  ionization  can  be  determined  by  considering 
only  the  ionization  of  sodium  by  atomic  oxygen. 

2.  For  altitiides  down  to  the  \  cinity  of  100,000  -  90,000  feet,  only  the 
forward  reactions  need  be  considered  Below  these  altitudes  the  sodium 
ionization  will  be  close  to  equilibrium. 

3.  Sodium  ionization  will  be  significant  with  respect  to  ciean  air  ionization 
depending  primarily  on  the  wall  mass  fraction  of  sodium  atoms.  A  mass 
fraction  of  10  ^  will  ensure  negligible  sodium  ionization. 

4.  Diffusion  works  to  reduce  the  peak  ionization  in  the  boundary  layer 
by  approximately  two  orders  of  magnitude,  whether  one  considers  clean  air 
or  contaminant  ionization. 

5.  The  relative  importance  of  clean  air  and  contaminant  ionization  is  a 
function  of  altitude  although  not  to  the  degree  that  it  is  a  function  of  wall 
mass  fraction  of  contaminant.  The  lov/er  the  altitude  the  more  dominant  the 
sodium  ionization.  Depending  on  the  wall  mass  fraction  of  sodium,  there 
will  be  an  altitude  above  which  the  clean  air  ionization  exceeds  the  sodium 
ionization. 


<  * 
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APPENDIX  A 


EQUILIBRIUM  IONIZATION 

The  non-equilibrium  analysis  assumes  that  the  ionization  of  sodium  is  far 
from  equilibrium.  Before  the  solution  for  any  particular  problem  is  accepted, 
the  validity  of  this  assumption  should  be  determined. 

The  equilibrium  between  an  atom  and  its  ions, 

A  -  A++  e"  {24) 

can  be  shown  to  be  described  by: 

1/2 

*  ”  P +P"K  i 


where  t  =  equilibrium  degree  of  ionization 

K  =  equilibrium  constant 
P 

P  =  pressure  (atmospheres) 

The  equilibrium  constant  is  described  by  the  equation  of  Saha®. 


logKp=^L+flogT-6.,  +  log  \SUg] 

where  E  =  ionization  potential  (electron  volts) 

T  =  temperature  (°K) 
g  =  statistical  weight 

The  last  term  can  be  considered  negligible^. 

The  electron  concentration  is  related  to  the  degree  of  ionization  by: 
N.-  -  P  « 


where  Ne  =  electron  concentration 

Cj^a  =  mass  fraction  of  sodium  atoms 
:=  atomic  weight  of  sodium 
p  =  mixture  density 
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K  rK  “i  l/2 

Using  P  =  p  R  T  and  £  =  =  p2  \  for 

_  p  J  _  I 

P  »  K  {  which  for  typical  peak  boundary  layer  temperatures  is 

It 

valid  for  altitudes  below  200,000  feet) 

Ne-  _  (Kpjl/2  CNa 
^UT  ~  RT 

The  local  mass  fraction  of  sodium  atoms,  C^,  is  directly  proportional 
to  the  wall  concentration.  Thus: 

Nv _  (kp)1/2 


RT 


where  the  right  side  is  a  function  of  temperature  only.  Thus  the  peak 

equilibrium  concentration  for  a  given  peak  temperature  is  cons  .ant  when 

Ne 

considered  in  the  form  — =■ -p* -  and  is  plotted  as  a  horizontal  line  in 

p  1  ^  Tr 

Figure  3.  ^  ^w 
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APPENDIX  B 


EFFECTS  OF  DIFFUSION 

It  is  of  interest  to  determine  how  much  effect  the  boundary  layer  diffusion 
process  has  in  reducing  peak  electron  concentration  levels.  The  prediction 
of  electron  concentrations  neglecting  diffusion  is  straightforward. 

The  conservation  of  mass  for  any  species  i  requires  that: 

div  {p  (q  r  qi)  ^  (30) 

v/here  q  =  mass  averaged  velocity 

q^  =  diffusional  velocity  of  species  i  with  respect  to  q. 

Considering  one -dimensional  flow  and  neglecting  diffusional  velocities,  this 
reduces  to: 


^  (P  U  Ci)  =  Wi  (31) 

which  integrates  to: 

d  =  Wi  ~~  ,  (32) 

considering  Wi  to  be  constant  (or  for  a  specific  value  of  q  only) 

If  we  consider  only  what  was  found  to  be  the  most  significant  reaction, 
namely  the  ionization  of  sodium  by  collision  with  atomic  oxygen,  then: 


w 


+  =  p^  C 


(cO)‘ 


and 


Na 


r  i  _  P  x 
CNa +  “  U 


Na  Mq^ 


(33) 


IF? 


k3  C 


Na 


ISO) 

[Moj 


(34) 


v/here  C^T  +  =  mass  fraction  of  sodium  ions 
Na 


U  sing 


N  =  N 
e  1NNa 


+  =  p  CNa+ 


M»t  + 
Na 


(35) 
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Equation  (34)  becomes: 


£  -  9.  (g| 

where  Ng“  -  electron  density  (electrons/cm^) 
x  =  distance  along  surface  (ft) 

R  =  gas  constant 

=  53.  3  £t  '  lbf 


lbm  -  OR 


Substituting,  and  converting  x  from  flat  plate  to  conical  surface  distances 

by  dividing  by  three,  ^ —  =  0.  35  (10)36  (CNa)  (Cq)2 

P  ^ 

It  should  be  noted  in  using  the  above  equation  that  Cn  varies  linearly  with 

Px  7  u 

•y — .  Levinsky  and  Fernandez  show  that  neglecting  diffusion 


„  _  Px  (3,58)  ( 1 0) 

"  U 

where  T  is  in  degrees  Kelvin. 


59,400 


Substituting  for  Cq,  Equation  (37)  becomes: 

Ne  _  „  „  ,,„*69  I  Px  |3  CNa.  -  212  600 


=  4.7  (10)' 


I  Px  |  CNj 
[U  I  T7 


where  T  is  in  degrees  Rankine. 

Levinsky  and  Fernandez  present  a  similar  result  for  the  clean  air  electron 

concentration,  namely 

.  .  !  Px \ 3 


^e  _  0.47  (10)48 
?  (T/ 1000) 


-  185,  000 


where  T  is  in  degrees  Kelvin. 
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(  '  ) 


APPENDIX  C 


(  •  ) 


RELATION  BETWEEN  WALL  BLOWING  AND  WALL  MASS  FRACTK'N 

If  it  is  assumed  that  the  velocity  of  air  within  and  up  to  the  surface  is 
zero* 


then  at  the  wall 
-  1 


.1 


<pv>w  =  t^  :;“37iw  («> 

where  C  here  refers  to  the  mass  fraction  of  total  material  injected  into  the 
boundary  layer. 

Transforming  to  {*1,  £)  where: 


Ue  ro 


(x)  ,y 

"1  =  — — ’TTT~  /  p  dy 
(2  i)l/Z  0 

£  =  Pe  ue  P-e  /*  ro  ds 

6 

and  introducing  the  stream  function 

*  (>i.5)  =  (2  5) 1/2  f(n) 


where  p  v  rQ  = 


,.x 


Then  Equation  (41)  becomes 


f(o)  = 


1 


1  *  SCi^ 


j 


(42) 

(43) 

(44) 

(45) 

(46) 


rj  ! 

'  !  W 

If  C  is  taken  to  be  non-reacting  so  that  its  normalized  solution  is  represented 


(47) 


by  Equation  (6),  then 

l  1 

j 

1 

„  i  /!  Sc 

jOC 

i  -  _  c 

os  ; 

^w  Lf  (o).i 

-  o 

w 

1  w 

69 

r  [£-(>1)1 Sc  a  >i 

o  u 

Taking  Sc  -  i  =  1,  Equation  (46)  becomes: 


f(o)  =  - 


Cw  f"(o) 


(48) 


w 


<  i  I 
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or  solving  for  Cw, 


r  -  *(o) 

-  {(of  -  FloT 

(49) 

Then 

r,  _  Q  f(0) 

~  f(o)  -  f"(o) 

(50) 

For  a  cone, 

£(o)  =  -V|  ■ 

yPe 

(51) 

where  if  the  solution  is  to  be  similar, 

(pv)w  '\/> c  =  constant 

(52) 
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NOMENCLATURE 


D 

f 

g 

H 

h 

h 

s 

Kw 

K 

P 

a 

M 

N. 

1 

P 

Pr 

Q 

R 


mass  fraction  of  species  i 
binary  diffusion  coefficient 
velocity  function  =  / yy-  d  q 

dimensionless  enthalpy  defined  by  equation  (17) 

viscous  dissipation  and  work  term  defined  by  equation  (20) 

thermal  enthalpy 


/. 


C  d  T 
Tre£  p 


total  enthalpy 


h  +  U2/2 


mass  fraction  of  sodium  atoms  at  the  wall 
equilibrium  constant 

forward  reaction  rate  constant  associated  with  reaction  i 
density  -  viscosity  ratio 

P^/Pe^ 

molecular  weight 

concentration  of  species  i  in  moles  per  unit  volume 

pressure 

Prandtl  number 

parts  per  million  of  sodium  contamination  of  injected  material 
divided  by  10^ 

gas  constant 
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S  normalized  atom  mass  fraction  defined  by  Equation  (2) 

Sc  Schmidt  number 

T  temperature 

U  velocity  parallel  to  surface 

v  velocity  normal  to  surface 

mass  rate  of  formation  of  species  i  per  unit  volume 

x  distance  along  surface 

€  equilibrium  degree  of  ionization 

p  mixture  density 

r|  transformed  independent  variable  (Equation  42) 

£  transformed  independent  variable  (Equation  43) 

SUBSCRIPTS 

e  edge  of  boundary  layer  conditions 

e  electron  species 

eq  equilibrium  conditions 

m  mean  or  average 

o  atomic  oxygen  species 

w  wall  conditions 

oo  free  stream  conditions 
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